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ABSTRACT 


In  this  article  the  dual  of  the  multi-facility  location  problem  with 

arbitrary  norms  is  developed.  The  formulation  allows  any  number  of  linear  f“r. 

constraints  in  the  primal.  It  is  shown  that  the  multipliers  associated  with  • TIS 

: VCD 

the  linear  equations  of  the  dual  are  the  optimal  facility  locations  of  the 
primal. 


AMS  (MOS)  Subject  Classitication : 90B10 

KeyWords:  Multi-facility,  location  problem,  dual,  arbitrary  norms, 
multipliers 


Work  Unit  Number  5 (Mathematical  Programming  and  Operations  Research) 


^Sponsored  in  part  by  the  United  States  Army  under  Contract  No.  DAAG29-7  5- 
C-0024.  This  research  was  supported  by  grants  from  the  Graduate  School  of 
the  University  of  Wisconsin:  use  of  the  University  of  Wisconsin  Computing 
Center  under  these  grants  was  made  possible  through  support,  in  part,  from  the 
National  Science  Foundation,  other  U.S.  government  agencies  and  the 
Wisconsin  Alumni  Research  Foundation. 

* 

Visiting  Professor  in  the  Department  of  Management  Sciences,  University 
of  Waterloo,  Waterloo,  Ontario  in  1976-77. 


ON  THE  DUAL  OF  THE  LINEARLY  CONSTRAINED  MULTI -FACILITY 
LOCATION  PROBLEM  WITH  ARBITRARY  NORMS1 
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We  consider  the  problem 


minimize 

subject  to  AjXj  £ b . 


ai 


) + l 


n-1 

i=l 


rn 

ij=i+l 


W2ijK2ij(xi 
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m and  n are  the  number  of  existing  and  new  facilities,  respectively,  and  d is 
the  dimension  of  the  facility  space.  is  the  non-negative  weight  and 

is  the  norm  on  R^,  to  be  used  between  existing  facility  i and  new 
facility  j for  1 £ i £ ®»  1 f j J n.  w2ij  *S  t*le  non_nc8at;^ve  weight  and 
^2ij  (" ) *s  t*le  norm  on  ***,  to  use<^  between  new  facilities  i and  j for 
1 £ i < j £ n.  a^  c R1*  Is  the  location  of  existing  facility  i for  1 £ i £ 
and  Xj  e R^  is  the  unknown  location  of  new  facility  j for  1 £ j £ n.  I linear 
constraints  on  the  locations  of  the  new  facilities  are  expressed  using  n l x d 
matrices  Aj  for  1 £ J £ n and  the  i-vector  b. 

This  article  generalizes  and  synthesizes  some  results  obtained  by  Love  [1] 
and  Planchart  and  Hurter  [4].  Using  elements  from  conjugate  function  theory 
as  cxpositcd  by  Witzgall  [5],  we  shall  find  the  dual  of  the  location  problem. 
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Let  superscript  c denote  the  conjugate  lunctlon  corresponding  to  a convex 


function,  and  superscript  o denote  the  polar  corresponding  to  a norm.  It  is 
T 


T T T 

(x. , x,,  * • • , x I and  A • (A. , A* , , ■ # , A ) • Then 
i-  £ n l r j n J 


convenient  to  let  x 

Ax  - A^Xj,  For  1 < i < u,  and  1 < J < n,  let  - (0,0, . . .0,1 ,0,0, . . . ,0) 

be  a d x nd  matrix  consisting  of  n d x d submatrices , each  of  which  is  a null- 

matrix,  except  the  J^h,  which  is  an  identity  matrix.  Then  x.  “ B,.,x.  For 

1 i J n J 11J 

1 i 1 K J ^.n,let  ®2ij  “ (0, 0, . . . ,0, 1,0,0,.  ..,0,-1, 0,0.  ..,0)  be  a dx  nd  matrix 

consisting  of  n d x d submatrices,  each  of  which  is  a null-matrix,  except  the 

1th,  which  is  an  identity  matrix  and  the  J.th,  which  is  a negative  Identity 

matrix.  With  this  notation,  we  can  write  down  a number  of  functions  occurring 

in  the  location  problem  and  find  the  conjugate  functions. 

For  1 < i < m,  and  1 < j < n,  let  f UJ  (x)  - <Blijx_a1) , so  we  have 


fliJ(Blijylij) 


3i  yliJ 


lf  Klij(ylij)  - Wlij 
otherwise 


For  1 < i < J < n let  f 2ij (x)  ” w2i jK 2i j *B2i jX* ’ S°  Ue  have 


f2iJ(BIijy2ij) 


if  K2ij(y2iJ)  - W2ij 


Let  f3(x) 


■>  otherwise 
if  Ax  < b 


• otherwise 


. Then  fjfcj)  “ sup{*  jX'.A^  <_  b)  - inf{bTy3  : ATy3=z3,y3>0) 


by  the  theory  of  linear  programming 

rn 


Finally,  let  f(x)  « ^ fUj(x)  + f (x)  + f3(x),  so 

the  location  problem  is:  minimize  f(x). 

For  the  dual  we  get: 


j 
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w* 


fc(0)  - a^y^  + inf(bTy3:ATy3  - ty  yy  >_  0): 

^1-l^j-l  Blij  yUJ  + ^i-1  ^J-i+1  B2iJ  y21J  + z3  “ °i  ~ wnj. 


1 < 1 < ».  l£  j <n;  K2ij(y2ij)  - U21J*  1£1<J£n}“ 

1"far.i  5-1  \ yUJ  + bTy3: 

^1-1  ylij  " ^1=1  y2ij  + ^i=j+iy21j  + A)y3  ” °»  1-J-n:  KliJ(yliJ)  - W11J’ 

1 l 1 £ m*  i i j < n;  K°ij^y2ij^  - w2ij*  1 l 1 * J £ n>  y3  < °). 

The  dual  of  the  location  problem  may  be  written  as 

. . r®  rn  T , T 

maximize  - ^ yuj  - b y3 

subject  to  ["„!  yUj  - l{l\  y2ij  + 5-J+1y2iJ  + A-V3  “ 0 for  1 < j < n, 

Kllj  ^yli  — Wlij  for  1 - 1 - m and  1 i 1 1 n> 

K21j(y2ij)  - W2ij  for  1 £ 1 < J 1 n» 

y3  > 0. 

This  dual  may  be  transformed  back  into  the  location  problem,  using  differentiable 
duality  theory  as  exposlted  by  Luenberger  [3].  Consider  the  problem  (P)  : 


maximize  p(y)  subject  to  g(y)  £ 0,  h(y)  = 0 (ycR  ) 

where  p(*)  is  a differentiable  concave  real-valued  function,  g(-)  is  a differ- 
entiable vector-valued  function  with  convex  component  functions  and  h(’)  is  an 
affine  function.  A dual  of  problem  (P)  is  this  problem  (D) : 


X 

minimize  d(x)  « max{p(y)  + x h(y):g(y)  £ 0). 

In  a similar  manner  to  the  development  by  l.uenbergcr  [3],  we  may  arrive  at  the 
following  facts: 
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1.  P(y)  £ d(x)  lor  all  y with  g(y)  <_  0,  h(y)  - 0 and  all  x. 

2.  Let  a constraint  qualification  hold  for  problem  (P) . If  y^  solves 

problem  (P)  with  multiplier  vector  corresponding  to  the  equality 

constraints,  then  x solves  proMem  (D)  and  p(y  ) - d(x  ). 

o o o 

These  facts  cannot  be  applied  directly  to  the  dual  of  the  location  problem, 
since  the  functions  in  the  nonlinear  constraints  are  norms  and  thus  not  dif- 
ferentiable. For  many  norms,  however,  a simple  transformation  will  yield  a 
differentiable  function.  Consider,  for  instance,  the  norm  K on  given  by 

K(x)  «*  t (C  ) where  C is  a non-singular  d x d matrix  and  l (z,,z_ z.)  = 

P * p 1 2 a 

(£i_l  The  Inequality  K(x)  £ w is  equivalent  to  (K(x))P  <_  wP,  in 

which  inequality  the  left-hand  side  is  a differentiable  function  for  p > 1. 

For  simplicity,  we  do  not  exhibit  such  transformations  explicitly  in  the  fol- 
lowing. 

Taking  the  dual  of  the  location  problem  as  (P) , we  obtain  for  problem  (D) : 
d(x)  = maxf-I^  XJ_X  *1  yUj  ~ + 

Ij-1  xj(Ii-i  ylij  ~ y2ij  + ^i-J+1  y2ji  + Ajy3): 

Klij<ylij)  - wlij  ’ 1 - 1 - m*  1 £ 3 5 n5  ^ij^aij5  £ w2ij » 1 - 1 < J - n»  y3^0j 

" C-i  £j-l  max{  (xj  " ai)Tylij  : Klij(ylij)  ±wlij}  + 

^i-1  ^j"i+l  max{(xi  “ xj3  y2ij  : K2ij(y2iJ*  -w2ij)  + 

"ax(^J-l  AJXJ  “ b)Ty3  : y3  1 0}* 

Using  the  properties  of  polars  of  norms  (5),  problem  (D)  is  thus  identical  to 
the  original  location  problem. 
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The  location  problem  can  be  Bolved  by  solving  lta  dual  by  some  standard 
nonlinear  programming  algorithm,  the  optimal  locations  being  the  optimal  mul- 
tiplier vectors  corresponding  to  the  equality  constraints.  Our  computational 
experience,  however,  seems  to  indicate  that,  from  the  viewpoint  of  computation 
time,  it  is  somewhat  more  efficient  to  solve  the  location  problem  directly, 
using  differentiable  approximations  to  the  norms  Involved  [2],  (Some  may  find 
the  dual  form  easier  to  program,  however,  due  to  the  simple  nature  of  the  de- 
rivatives of  the  objective  function  and  constraints.)  The  special  structure 
of  the  dual  can  be  exploited  to  solve  it  using  decomposition  as  suggested  by 
lx>ve  [1)  and  Planchart  and  Hurter  (4J.  Even  this  technique  for  solving  the 
location  problem  seems  less  efficient  than  either  solving  the  primal  form  or 
solving  the  dual  directly  with  a nonlinear  programming  routine. 

ACKNOWLEDGMENTS 

This  research  was  supported  by  the  Graduate  School  and  the  Graduate  School 
of  Business  of  the  University  of  Wisconsin,  Madison.  Use  of  the  University 
of  Wisconsin  Computing  Center  was  made  possible  through  support,  in  part,  from 
the  National  Science  Foundation,  other  United  States  government  agencies,  and 
the  Wisconsin  Alumni  Research  Foundation  through  the  University  of  Wisconsin 
Research  Committee. 

REFERENCES 

1.  R.  F.  LOVE,  "The  Dual  of  a Hyperbolic  Approximation  to  the  Generalized 

Constrained  Multi-Facility  Location  Problem  with  Ip  Distances,"  Manage- 
ment Sci.  21,  22-33  (1974). 

2.  AND  J.  C.  MORRIS,  "Solving  Constrained  Multi-Facility  Location  Prob- 

lems Involving  ip  Distances  Using  Convex  Programming,"  Opns.  Res.  23, 
581-587  (1975). 


-5- 


3.  DAVID  C.  LUENBERCER,  Introduction  to  Linear  and  Nonlinear  Programming, 

Addison-Wesley , Reading,  Massachusetts,  1973,  3X2-316. 

4.  A.  PLANCHART  AND  A.  P.  HURTER,  JR.,  "An  Efficient  Algorithm  for  the  Solu- 

tion of  the  Weber  Problem  with  Mixed  Norms,"  SIAM  J.  Control  13. 
650-665  (1975). 

5.  C.  WITZGALL,  "Optimal  Location  of  a Central  Facility — Mathematical  Models 

and  Concepts,"  National  Bureau  of  Standards  Report  No.  8388  (1964). 


-6- 


TITLE  (*nd  Subtitle) 


j (j  numbi :n 


ICCUNITV  u.  MSiriCAIIf'M  O’  THIS  >>»<.!  f’etm  linletej)  

| REPORT  nOCUWCin  ATIOM  PAGE nr.roNi:  cc^iin.rTiNo  i ohm  ) 

I!  nfPom  M'u'ir*  ' — . ji.  oovt  Acce»»io^fTO.  » i cat al4,u  HUHttn  J 

j_72j ^ ch  n I *c <a  I \ i 

4 . TITLE  (v<1  lufc'iif*'  | ]|  COVLRtO 

\£j  £oN  THE  DUAL  OF  THE  LINEARLY  CONSTRAINED  I I nummary /CpW~^  no  specific 

j MULTI  -FACILITY  LOCATION  PROBLEM  WITH  r 

ARBITRARY  NORMS.  I 


6 PIIIFORMINO  ORC 


" rt-muu  i.uvu(tu 

t.V  no  specific 

mod  ^ 

REPORT  NUMBER 


[7.  AUTHORS#) 


B.  CONTRACT  OR  GRANT 


Henrikyjuel  off B Robert^Love 


DAAG29-7  5-C-0/24 


8 PERFORMING  ORGANIZATION  NAME  UNO  AOOnESS  '0.  ^ FI  AM  FX  EM  EN  T.  FMIO  J E C T , T ASK  j 

_ _ , __  . ..  r AKtA  O WOHK  UNIT  NUMULRj  I 

Mathematics  Research  Center,  University  of  y 5 (Mathematical  Programming  ! 
610  Walnut  Street  Wisconsin  and  Operations  Research)  • 

Madison.  Wisconsin  53706 , 

II.  CONTROLLING  OFFICe  NAME  and  ADDRESS  S \I2.  ^RCPOHT  HATF  —f  ' 


tu.  S.  Army  Research  Office 
. P.O.  Box  12211 

i Research  Triangle  Park.  North  Carolina  27709 


y/£js 


MONITORING  \GF.NCY  MAME  A ADORESS(/f  dlliermnt  from  Controlling  Olllce) j '5-  SECURITY  CLASS,  (of  thle  report)  [ 


0J)% 

— ^ — /I 


UNCLASSIFIED 


ls«.  OECI.  ASSI  FI  CATION/ OOWN  OR  AGING 
SCHEDULE 


UT  DISTRIBUTION  STATEMENT  (ol  title  Report, 


W}\RC~' T$R-d73/} 


f Approved  for  public  release;  dlslH5ution  unlimited. 


I?  DISTRIBUTION  STATEMENT  (o t the  abetted  entered  In  Block  30,  It  dllletent  Irom  Report) 


l l»  SUPPL EMEN T ARY  NOTE 


• t.  KEY  WORDS  (Continue  on  reeeree  elde  II  neceeeery  end  Identity  by  block  number) 

Multi-facility 
Location  problem 
Dual 

Arbitrary  norms 

•'Multipliers 

JO.  ABSTRACT  (Continue  an  revet  ee  el  Jo  II  neceeeery  mnd  Identity  by  block  murtber) 

• In  this  article  the  dual  of  the  multi-facility  location  problem  with 

arbitrary  norms  is  developed.  The  formulation  allows  any  number  of  linear 
constraints  in  the  primal.  It  is  shown  that  the  multipliers  associated  with 
linear  equations  of  the  dual  are  the  optimal  facility  locations  of  the  primal. 


DD  1473  coition  or  i nov  «s  is  obsolete 


UNCLASSIFIED 


SrClIfilTV  O AMlFir  ATION  OF  TmH  pa  Of  (When  Hate  entered 


